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In this work partition of spaces algorithm is proposed to solve optimal 
reactive power problem. In this algorithm, for finding the finest outcome 
based on the concentration of elevated quality and capable points in specific 
area is considered. State space area are identified and divided into subspaces 
iteratively and search has been made more comprehensively. Performance of 
the proposed partition of spaces algorithm is evaluated in standard IEEE 
118,300 bus systems and simulated outcome gives better results. Real power 
loss has been considerably reduced. 
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1. INTRODUCTION 

Optimal reactive power problem has been key problem in power system, since it plays major role in 
secure and economic operation of the power system. Many conventional methods [1-6] have been applied for 
solving optimal reactive power problem. But many drawbacks have been found in the conventional methods 
and mainly difficulty in handling the inequality constraints. Last two decades many evolutionary algorithms 
[7-18] continuously applied to solve the problem. In this paper, partition of spaces algorithm is proposed to 
solve the reactive power problem. In this approach, elevated quality and capable points of the area is taken. 
State space are identified and divided into subspaces iteratively and search has been made more 
comprehensively. Performance of the proposed algorithm is evaluated in standard IEEE 118,300 bus systems 
simulation results shows the better performance of the proposed algorithm in reduction of real power loss. 


2. PROBLEM FORMULATION 

The key objective of the reactive power problem is to minimize the system real power loss 
and given, 

Ploss= £ k=l Sk(V?+V?-2Vi Vj cos6ij ) 
k=(ij) 1J 

Voltage deviation magnitudes (VD) is stated as follows: 
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Minimize VD = — 1.01 

(2) 

Load flow equality constraints: 


' G m cos Oil 

Pci - Pdi - v i Vj +Bj . sin 0.. =0u= 1,2 ...., nb 

(3) 

' Gjj sin 0;; 

Qgi - Qdi - V i Vj -t-Bjj cos 0jj = °’ 1 = 1,2 - ■ ’ nb 

(4) 

Inequality constraints are: 


V G T< V Gi <V£ ax ,iEng 

(5) 

V L r < V Li < V L ™ ax , i £ nl 

(6) 

Q£ in < Qci<Q£ ax ,ienc 

(7) 

Qm" < Qgi < QSfM G ng 

(8) 

T.min < T . < T max ; j £ nt 

(9) 

s min < gmax^ j £ n j 

GO) 


3. PARTITION OF SPACES ALGORITHM 

In this algorithm, for finding the optimal solution based on the concentration of elevated quality and 
capable points in specific area is considered. With equal sizes the state space are alienated into some 
subspaces. In the state space uniformly, points are generated arbitrary mode and target function value is 
calculated. Promising points are chosen with each subspace is determined. Always the chances of finding 
the optimal solution are higher when promising points are considered as the promising subspaces. The details 
of the proposed algorithm are as follows: 

a. Initially whole state space is measured as the capable area. 

b. Subspaces are created from the state space. Until finally a grid containing gi-L X gi 2 X ...gi n subsections 
in d th dimension (1 < d < n) and is divided into as many as g d equivalent subintervals. Then, Size of 
each subinterval, 

in dth dimension = (11) 

gid 

c. Preliminary population are created arbitrarily and that population is considered as the existing 
population. 

d. For every point of the existing population of function “f ’ is computed. 

e. Q1 - % of points; existing population are regarded as capable points, with the lowest values of function 
f, 

f. In each of the subspaces the numbers of capable points are determined and indicate the degree of 
capable points in the subspace. 

Capable ranks = Number of capable points in “s” (12) 

g. Q2 subspaces are searched more accurately and comprehensively. 

1) Extra specific search in the capable subspaces: smaller subspaces are made from the capable 
subspaces. 

2) To search subspaces expansively the number of points generated in each subspace is a linear 
function and found by [19]: 


number os points generated in s 


( capable rank s 
Steal! sub spaces capable rank k xpopulation size 


) 


( 13 ) 
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h. For each point in the new population the value of function f is calculated. 

i. Based on the truncation selection new population are replace the current population. 

j. Q3 percent of the points are arbitrarily selected in new population and customized by adding a Gaussian 
noise to them it’s alike to mutation operation in genetic algorithm. 

k. Evaluation of the stop condition. 

l. The output of algorithm is the most excellent solution generated so far. Trade-off between exploration 
and exploitation is controlled by the values of parameters Ql, Q2, and Q3. 

Partition of spaces algorithm for solving reactive power problem: 

Exploration Space, f 
{Inputs: 

With boundaries the n-dimensional space - State Space 
Function target denoted by "f" 

Most excellent solution found by the algorithm is the output 

Pop Size, Ql, Q2, Q3, gd for 1 < d ^ n are initialized 

When PopO = the initial population; i = 0;points are generated. 

Whole exploration space as the skilled area at the commencement; competent Subspaces= State 
Space; 1 < d ^ n for any dimension d 
State Space is Partition into gid parts; 

When end condition not satisfied 

{ 

Values of Good Points; Quality [,... Pop Size] = the quality of all the points in Popi 
Ql points with the uppermost superiority = [1... Pop Size] 

Number of superior points in each subspace Count = [1, . . . Number of Subspaces] 

count [1...number os sub spaces] , , , 

-= capable rank II... number os sub spaces] 

pop size 

Q2 percent of the subspaces with highest capable rank = Capable Sub spaces 

Dimension d : 1 < d ^ n for any proficient subspaces "s" ; Promising subspace "s" is 

partitioned into gd parts; 

Fresh generated subspaces = popj[l... pop size] 

Most excellent points of 
Popi = pop i+1 [l... pop size]; 

Gaussian noise is added to Q3 arbitrarily and chosen Popi+1 ; i = i +1; 

} 

Superior solution found so far will be the out put 
Revert to Solution 
} 


4. SIMULATION RESULTS 

Performance of partition of spaces algorithm evaluated in standard IEEE 118-bus test system [20], 
Limitations on reactive power source are listed in Table 1. Table 2 clearly shows better performance of 
the proposed algorithm. In IEEE 300 bus system, partition of spaces algorithm is verified [20]. Table 3 shows 
the values of comparison. Figure 1 gives real power loss comparison for IEEE 118 bus test system and 
Figure 2 shows the active power loss comparison for IEEE 300 bus test system. 


Table 1. Reactive power sources limits 


BUS NO 

5 

34 

37 

44 

45 

46 

48 

74 

79 

82 

83 

105 

107 

110 

QC 

MAX 

0.00 

14.0 

0.00 

10.0 

10.0 

10.0 

15.0 

12.00 

20.0 

20.0 

10.0 

20.0 

6.00 

6.00 

QC 

MIN 

-40. 

0.00 

-25. 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


Table 2. Real power loss comparison results 


Active power loss 
(p.u) 

Method- BBO 
[21] 

Method- ILSBBO/strategyl 
[21] 

Method -ILSBBO/Strategy2 
[21] 

Projected algorithm - 
Partition of spaces algorithm 

Minimum 

128.770 

126.980 

124.780 

116.640 

Maximum 

132.640 

137.340 

132.390 

120.980 

Average value 

130.210 

130.370 

129.220 

117.560 


Table 3 Comparison of loss 


Parameter 

Method CSA 

Method EGA 

Method EEA 

Partition 


[22] 

[23] 

[23] 

of spaces 

PLOSS (MW) 

635.8942 

646.2998 

650.6027 

629.5431 
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Figure 1. Real power loss comparison 


Figure 2. Active power loss comparison 


5. CONCLUSION 

Partition of spaces algorithm has been efficiently applied for solving reactive power problem. In this 
approach, elevated quality and capable points of the area is taken. State space are identified and divided into 
subspaces iteratively and search has been made more comprehensively. Efficiency of the proposed algorithm 
is evaluated in standard IEEE 118,300 bus systems and simulated outcome gives better results when 
compared to other reported standard algorithms. 
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